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ABSTRACT 

We discuss the 6 Li abundance evolution within a hierarchical model of Galaxy for- 
mation which correctly reproduces the [Fe/H] distribution of metal-poor halo stars. 
Contrary to previous findings, we find that neither the level ( 6 Li/H= 6 x 10 -12 ) 
nor the flatness of the 6 Li distribution with [Fe/H] can be reproduced under the 
most favourable conditions by any model in which 6 Li production is tied to a (data- 
constrained) Galactic star formation rate via cosmic ray spallation. Thus, the origin of 
the plateau might be due to some other early mechanism unrelated to star formation. 
Key words: nuclear reactions, nucleosynthesis, abundances - stars: abundances - 
stars: formation - cosmic rays - galaxies: evolution - cosmology: theory. 



1 INTRODUCTION 

The relative abundance of light elements synthesized during the 
big bang nucleosynthesis (BBN) is a function of a single parame- 
ter, r\, namely the baryon-to-photon ratio. Given the WMAP con- 
straint r\ = (6.8 ± 0.21) X 10~ 1 , the light nuclei abundances ca n 
be precisely predicted by BBN dSpergelfl2007l : lYao fe et~al1l2006l ). 
Despite a general agreement with the observed abundances of 
light elements, discrepancies arise concerning Li abundance. Ob- 
servationally, the primordial abundance of lithium isotopes ( 7 Li 
and 6 Li), is measured in the atmospheres of Galactic metal-poor 
halo stars (MPHS). 

bince the first detection bvlSpite fe Spitel fll982T), later con- 
firmed by subsequent works |Spite et alHl98*3 1 Ryan et al .11 19991 ; 
lAsplund et alj|2006t lBonifacicll2007ft a 7 Li/H = (1 - 2) x 10" 10 
abundance was deduced, independent of stellar [Fe/H]. The pres- 
ence of such a 7 Li plateau supports the idea that 7 Li is a pri- 
mary element, synthesized by BBN. The measured value, how- 
ever, results of a factor 2 — 4 lower than the exp ected from 
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, ^ 1Q- 10 JCuoco et alj[2004), or 7 Li/H = 4 . 15±° - f 5 X 
_ ~ 10 lICoc et al.ll2004T) . Recently, Ipinsonneault et al.l <2002ft and 
iKorn et^lTl2*00a ) found that mixing and diffusion processes dur- 
ing stellar evolution could reduce the 7 Li abundance in stellar 
atmospheres by about 0.2 dex, thus partially releasing the ten- 
sion. 

A more serious problem arose with 6 Li, for which the BBN 
predicts a value of ( e Li/H)BBN ~ 10~ 14 . Owing to the small 
difference in mass between 6 Li and 7 Li, lines from these two iso- 
topes blend easily. The detection of 6 Li then results quite dif- 
ficult since the predominance of 7 Li. Recently, high-resolution 
spectro scopic observations measured the 6 Li abundance in 24 
MPHS iAspIund et alj|2006h , revealing the presence of a plateau 
6 Li/H= 6 X 10" 12 for -3 £ [Fe/H] < - 1. A primordial origin of 
6 Li seems favoured by the presence of the plateau; however, the 
high 6 Li value observed cannot be reconciled with this hypothesis. 

The solutions invoked to overcome the problem were: 
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a modification of BBN mode l s llKawa saki et all 
Jcdam zlk et al.l 1200! | iPospelovl 120071 ; ICumberbatch et al.l 



2005 



2007 



Kusakabe et al.ll2007T) . (ii) the fusio n of 3 He accelerated by stella r 



flares with the atmospheric helium (Tatischcff & Thibaud 2007), 
(iii) a mechanism allowing for later production of 6 Li during 
Galaxy formation. The latter scenario involves the generation of 
cosmic rays (CRs). 6 Li, in fact, can be synthesized by fusion re- 



actions (a + a 



3 Li) when high-energy CR particles collide 



with the ambient gas. Energetic CRs can either be accelerated by 
shock wa ves produced during cosmological structure formation 
proce sses llMiniati et al.l2000l ; TSuzuki fc Inouell2002HKeshet et al.l 
l2003t) or, by strong supernova (S N) shocks along the b uild-up of 
the Galaxy. In their recen t work iRol linde e t alj | |200(J) used the 
supernova rate (SNR) by Dai gne et al.l (I2006T) to compute the 
production of 6 Li in the intergalactic medium (IGM). Assuming 
that all MPHS form at z ~ 3, and from a gas with the same 
IGM composition, they obtained the observed 6 Li value. Despite 
the apparent success of the model, these assumptions are very 
idealized and require a closer inspection. We revisit the problem 
using a more realisti c and data-cons t rained approach, based on 
the recent model by ISalvadori et al.l J2007t) (SSF07), which fol- 
lows the hierarchical build-up of the Galaxy and reproduces the 
metallicity distribution of MPHS. 



2 BUILDING THE MILKY WAY 

The code GAlaxy Merger Tree fc Evolution (gamet e) described 
in SSF07 (updated version in I Salvadori et al]|2008f) follows the 
star formation (SF)/chemical history of the MW along its merger 
tree, finally matching all its observed properties. 

The code reconstructs the hierarchical merger history of 
the MW using a Monte Carlo algorithm based on the ex- 
tended Press & Schechter theory llPress fc Schech tedll97 4) and 



adopting a binary sc heme with accretion mass (|Cole et al.l 2000l : 
IVolonteri et alj [2003). Looking back in time at any time-step a 
halo can either lose part of its mass (corresponding to a cumu- 
lative fragmentation into haloes below the resolution limit M IeB ) 
or lose mass and fragment into two progenitors. The mass below 
M res accounts for the Galactic Medium (GM) which represents 
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Figure 1. Upper panel: Comoving SFR density evolution for 
Pop III (solid line) and Pop II/I stars (dashed line). The curves 
are obtained after averaging over 100 realizations of the merger 
tree; shaded areas denote ±1<t dispersion regions around the 
mean. Poin ts represent the low-redshift measurements of the cos- 
mic SFR bv lHopkinsI ||2004| ) . Lower panel: Corresponding GM iron 
(solid line) and oxygen (dashed line) abundance evolution. The 
po int is the measured [ O/H] abundance in high- velocity clouds 
bv lGangulv et alj <2005l) . 



the mass reservoir into which haloes are embedded. During the 
evolution, progenitor haloes accrete gas from the GM and viri- 
alize out of it. We assume that feedback suppresses SF in mini- 
haloes and that only Lya cooling haloes (T v j r > 10 4 K) con- 
tribute stars and metals to the Galaxy. This motivates the choice 
of a resolution mass A/ rcs = M 4 (z)/10 = M(T vil = 10 4 K , z)/10 
where M±(z) is the mass corresponding to a virial temperature 
T v ir = 10 4 K at redshift z. At the highest redshift of the simula- 
tion, 2 w 20, the gas present in virialized haloes, as in the GM, 
is assumed to be of primordial composition. The SF rate (SFR) 
is taken to be proportional to the mass of gas. Following the 
critical met a llicity scenario llBromm et alj|200ll; ISchneider et al.l 
|2002|, 120031 ; lOmukai et ail |2005| ; ISchneider et al.l |2006|) we ~as- 
sume that low-mass (Pop II/I) SF occurs when the metallicity 
Z CT > 1O~ 5±1 Z0 according to a Larson initial mass function 
with a characteristic mass m* = O.35Af0. At lower Z massive 
Pop III stars form with a characteristic mass mp op in = 200A/q , 
i.e. within th e pair-instability supern ova (SN 77 ) mass range of 
140-260M Q llHeger fe Wooslevll2002l) . The chemical evolution of 
both gas in proto-Galactic haloes (ISM) and in the GM, is com- 
puted by according t o a mechanical feedback prescription (see 
ISalvadori et al.l | |2008|) for details). Produced metals are instan- 
taneously and homogeneously mixed with the gas. 

The model free parameters are fixed to match the global 
properties of the MW and the Metallicity Distribution Function 
(MDF) of MPHS derived form the Hamburg-ESO Survey (Beers 
& Christlieb, private communication). In Fig. [T] (upper panel) the 
derived Galactic (comoving) SFR density is shown for Pop III and 
Pop II/I stars. Pop II/I stars dominate the SFR at any redshift. 



Following a burst of Pop III stars, in fact, the metallicity of the 
host halo raises to Z > Z CT : chemical feedback suppresses Pop III 
formation in self-enriched progenitors. Later on Pop III stars can 
only form in those haloes which virialize from the GM and so, 
when Zqm Z cr , their formation is totally quenched. The above 
results are in agreement with recent hy drodynamic simulation s 
implementing chemical feedback effects l lTornatore et all |2007|) . 
The earlier Pop III disappearance of our model (2 ~ 10) with 
respect to this study (2 ~ 4) is a consequence of the biased vol- 
ume we consider i.e. the MW environment. As the higher mean 
density accelerates SF/metal enrichment, PopIII stars disappear 
at earlier times; the SFR maximum value and shape, however, 
match closely the simulated ones. 

In Fig. [T] (lower panel) we show the corresponding evolution 
of the GM iron and oxygen abundance. As SSF07 have shown 
that the majority of present-day iron-poor stars ([Fe/H]< —2.5) 
formed in haloes accreting GM gas which was Fe-enhanced by 
previous SN explosions, the initial [Fe/H] abundance within a halo 
is set by the corresponding GM Fe-abundance at the virialization 
redshift. 



3 LITHIUM PRODUCTION 

To describe the production of 6 Li for a continuous s ource of CRs 
we generalize the classical work of iMontmerlel |l977), who devel- 
oped a formalism to follow the propagation of an homogeneous 
CR population in an expanding universe, assuming that CRs have 
been instantaneously produced at some redshift. 

Since the primary CRs are assumed to be produced by SNc, 
the physical source function Q(E, 2) is described by a power law 
in momentum: 



Q(E,z) = C(z)^ (GeV/nrW 



with P ■ 



HE) 



- v/c and 
E-\ 



So 



[E(£ + 2£ )] (7+1)/2 



(GeV/n) -1 cm" 



(1) 



(2) 



where 7 is the injection spectral index and Eo = 939 MeV and 
E are, respectively, the rest-mass energy and the kinetic en- 
ergy per nucleon. The functional form of the injection spectrum 
(E) is inferred from the th eory of collisionless shock acceleration 
Blandford fc Eichlenll987l) and the 7 value is the one typically 
associated to the case of strong shock. We note however that the 
results are only very weakly dependent on the spectral slope. Fi- 
nally, C(z) is a redshift-dependent normalization; its value is fixed 
at each redshift by normalizing Q(E, 2) to the total kinetic energy 
transferred to CRs by SN explosions: 



£sn(z) 



with 



EQ{E,z)dE 



(3) 



(4) 



£ S nW = e(l + 2) 3 [E n SNR n (2) + , 

where En = 1.2 X 10 51 erg and E^ = 2.7 X 10 52 erg are, respec- 
tively, the average explosion energies for a Type II SN (SNjj) and 
a SN-yy; e = 0.15 is the fraction of the total energy not emitted 
in neutrinos transferred to CRs by a single SN, assumed to be 
the same for the two stellar populations; SNRn (SNR-yy) is the 
SNjj (SNyy) explosion comoving rate, simply proportional to the 
Pop II/I (Pop III) SFR. The efficiency parameter is inferred by 
shock acceleration theory and confirmed by re cent observations 
of SN remnants in our Galaxy <Tatischeffll2008ft . 



We now need to specify the energy limits E n 



of 



the CR spectrum produced by SN shock waves (eq. |3j. We 
fix £max — 10^ GeV, following the theoretical estimate by 
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lLagage fc Cesarskvl J 1983f> ■ Due to the rapid decrease of <t>(E) the 
choice of -E max does not affect the result of the integration and 
hence the derived C(z) value. On the contrary C(z) strongly de- 
pends on the choice of i? m i n : the higher -E m i n , the higher is C(z). 



Since observations cannot set tight constraints on E n 



due to 



solar magnetosphere modulation of low-E CRs, we consider it as 
a free parameter of the model. 

Once the spectral shape of Q(E, z) is fixed, we should in prin- 
ciple take in account the sub sequent propagation of CRs both in 
the ISM and GM. Following IRollinde et ail [|2006h . we make the 
hypothesis that primary CRs escape from parent galaxies on a 
timescale short enough to be considered as immediately injected 
in the GM without energy losses. At high redshift in fact: (i) 
structures are smaller and less dense (Zhao et al. 2003) imply- 
ing higher diffusion efficiencies (Jubelgas et al. 2006); (ii) the 
magnetic field is weaker and so it can hardly confine CRs into 
structures. Note also that, besides diffusive propagation of CRs, 
superbubbles and/or galactic winds could directly eject CRs into 
the GM. 

Under this hypothesis the density evolution of primary CRs 
only depends on energy losses suffered in the GM. The nuclei lose 
energy mainly via two processes, ionization and Hubble expan- 
sion, and they are destroyed by inelastic scattering off GM targets 
(mainly protons). 

We can follow the evolutio n of ct-particles (primary CRs) 
through the transport equation llMontmerlelll977f) 



at oE Tr, 



(5) 



where Ni h is the ratio between the (physical) number den- 
sity of species i and GM protons, nn(z) = riH,o(l + ^) 3 ; 
Q^n(E,z) = Q(E, z)/nn(z) is the normalized physical source 
function, b = (dE/dt ) is the total energy loss rate adopted from 
iRollinde et al. I ll2006h, Tp is the dest ruction term as in the ana- 
lvtic fit by iHeinbach fc Simonl jl995h : finally, K ap = 0.08 is the 
cosmological abundance by number of a-particles with respect to 
protons. 

We consider 6 Li as entirely secondary, i.e. purely produced 
by fusion of GM He-nuclei by primary a-particles. The physical 
source function for 6 Li is given by: 



hu(E,z) 



Li (E,E')n Ke (z)$ a (E',z)dE' 



(6) 



where E' and E are respectively the kinetic energies per nu- 
cleoli of the incident particle and of the produced 6 Li nuclei, 
and § a (E',z) = f3(E')N a (E' , z) the incident a-particle flux. 
Making the approxima tion a a a^6Li(E, E') = ai(E)8{E — E'/4) 
llMeneguzzi et al]|l97ll) and defining Q6Li,H = QeLi/ n Hj the eq. 
lO becomes 



Q6Li, H (E, z) = a t (E)K ap n H (z)*a,H (4£, z) 



(•) 



where the cross section a - ; (E) is given by the analytic fit of lMercerl 
(l200lh : 



odE) ~ 66exp ( — ) 

K ' \ 4MeV/ 



nib 



(8) 



We can now write a very simple equation describing the evo- 
lution of 6 Li: 



£W6Li,H 

at 



hu,s(E,z) 



(9) 



in this case, in fact, destruction and energy losses are negligible 
since their time scales are very lo ng with respect to the production 
time scale llRollinde et al.ll2005ri . 

The solution of the coupled eqs. J5jl-(|9jl gives 6 Li/H at any 
given redshift z. 



4 RESULTS 

The system of equations introduced in the previous Sec. are 
solved numerically using a Crank-Nicholson implicit numerical 
scheme jPressll2002T) . Because of its stability and robustness im- 
plicit schemes are u sed to solve transport equat ions in most CRs 
diffusion problems l lStrong fc M oskalcnko 1998). 

We test the accuracy of our code by studying a simplified 
case in which an analytic solution can be derived and compared 
with numerical results. To this aim we assume that: (i) both en- 
ergy losses and destruction of primary CRs in the GM can be 
neglected; (ii) the physical energy density injected by SNe is con- 
stant, £sn ~ 7.4 X 10 -27 GeV cm -3 s _1 , in the redshift range 
Z. > 3. It is worth noting that the above hypothesis conspire to 
give an upper limit to the exact solution, thus providing an es- 
timate of the maximum achievable 6 Li abundance. Under these 
approximations, the source spectrum defined in eq. Jl} becomes: 



(10) 



(11) 



Q(E, z) = 6.4 • 10" 29 4^ (GeV/n)" 1 cm" 3 s" 1 
/3(E) 

and eqs. J5]l-(|9]l can be solved. We find 
N a , H (z) = 39.6(1 + z)S/ 2 
and 

A r 6Li, H (^) = 8.2 X IQ-^il + z)- 3 (12) 

From Fig. [2] we conclude that the analytical solution for the 
GM 6 Li abundance (eq. I12H is perfectly matched by the numer- 
ical one. Also shown are the numerical solutions obtained by 
relaxing first the hypothesis (i) and then (i) + (ii). Not unex- 
pectedly, the inclusion of energy losses and destruction term into 
eq. {(5} affects only slightly the result, as the typical time-scales 
of such processes are longer than the 6 Li production one. 

A realistic injection energy, on the contrary, has a strong 
impact on the predicted shape and amplitude of the 6 Li evolution. 
In fact, the SFR, and consequently £gN, is an increasing function 
of time in the analyzed redshift range (Fig. 1 upper panel). The 
maximum we can obtain by using the SNR derived from the 
curve in Fig. 1, a r ealistic energy transf er efficiency e = 0.15, and 
_B min = 10" 5 GeV llRollinde et al.l2006h . is £S^ X ~ 8.6xl0" 28 < 
7.4 X 10" 27 GeV cm -3 s -1 . Note that the °Li/H abundance at 
Z = 3 results more than 1 order of magnitude smaller than the 
value of the simplified case. In the following, we will refer to this 
physical model as our fiducial model. 

We now use the [Fe/H] predicted by GAMETE (Fig. 1, lower 
panel) to convert redshift into [Fe/H] values and derive the GM 
6 Li vs [Fe/H]. According to our semi-analytical model for the 
build-up of the MW, in fact, the GM elemental abundances reflect 
those of MPHS, which are predicted to form out of new virializing 
haloes accreting gas from the GM. This implies that the observed 
MPHS formed continuously within the redshift range 3 < z 10. 
From Fig. [3] we see that our fiducial model yields log 6 Li/H= 
— 13.5, i.e. about three orders of magnitude below the data. 

This discrepancy cannot be cured by simply boosting the 
free parameters to their maximum allowed values. This is also 
illustrated in the same Figure, where for the upper curve we as- 
sume e = 1, -E m in = 10 MeV/n0 and for the SFR the maximum 
value allowed by GAMETE within l-a dispersion. Although the 
discrepancy between observations and model results is less promi- 
nent in this case, we are still unable to fit the data, in particular 



1 Thi s solution represents an upper limit for the IRollinde et ahl 
ll2006h model, as inferred from their Fig. 2. 

2 This value is exceptionally high and corresponds to the energy 
at which the 6 Li production is most efficient. Thus the 6 Li pro- 
duction will be drastically reduced by increasing E m - ln above this 
value. 
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Figure 2. Redshift evolution of GM 6 Li/H abundance for the 
analytical (green short-dashed line) and numerical solution (over- 
lapped) of a simplified model with no energy losses and de- 
struction, and £g N = 7.4 X 10~ 27 GeV cm -3 s _1 for z > 3, 
the same model including energy loss/destruction (blue long 
dashed line), the fiducial model with realistic SNR, e = 0.15 and 
E min = 10" 5 MeV (red solid line). 

at [Fe/H] = -3 (i.e. at higher redshifts) only Log 6 Li/H= -12.6 
has had time to be produced, failing short by 30 times. 

In addition the fiat data distribution canno t be recovered. 
It is w orth noting that, as also pointed out by lAsplund et al.l 
(2006) 6 Li may be depleted in stars, mainly during the pre-main 
sequence phase. If this is the case, the 6 Li abundance observed 
in stars would not be representative of the gas from which they 
have formed. Taking into account this effect the inferred 6 Li abun- 
dances become metallicity dependent, i.e. the flatness is lost. Be- 
cause of depletion however, the derived 6 Li values would be higher 
for all [Fe/H], making the discrepancy between our results and 
observations even larger. 

We finally note, as already claimed by Rollindc et al.l J2006T) ■ 
that the production of 7 Li through this mechanism is comparable 
with that of 6 Li , being the production cross sections of the two 
isotopes very similar. No overproduction of 7 Li is then expected 
with respect of the BBN-based value. 



5 DISCUSSION 

We have pointed out that both the level and flatness of the 6 Li 
distribution cannot be explained by CR spallation if these par- 
ticles have been accelerated by S N shocks inside MW building 
blocks. Although previous claims llRollinde et al.l l2006') of a pos- 
sible solutiorjfl invoking the production of 6 Li in an early burst of 
PopIII stars have been put forward, such scenario is at odd with 
both the global properties of the MW and its halo MPHS. 

Our model, which follows in detail the hierarchical build- 
up of the MW and reproduces correctly the MDF of the MPHS, 

3 Note that their eq. 18 contains an extra dz/dt term 
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Figure 3. Redshift evolution of 6 Li/H vs [Fe/H] for the fiducial 
model (e = 0.15, E mln = 10 GeV/n, dashed line) and for the 
maximal model (e = 1, E m i n = 10 MeV/n, solid line). Shaded 
areas denote ±1ct dispersion regions around the mean. 

predicts a monotonic increase of 6 Li abundance with time, and 
hence with [Fe/H], Moreover, our fiducial model falls short of 
three orders of magnitude in explaining the data; such discrep- 
ancy cannot be cured by allowing the free parameters (i£ m i n ,e) 
to take their maximum (physically unlikely) values. Apparently, 
a flat 6 Li distribution appears inconsistent with any (realistic) 
model for which CR acceleration energy is tapped from SNe: if 
so, 6 Li is continuously produced and destruction mechanisms are 
too inefficient to prevent its abundance to steadily increase along 
with [Fe/H]. 

Clearly, the actual picture could be more complex: for ex- 
ample, if the diffusion coefficient in the ISM of the progenitor 
galaxies is small enough, 6 Li could be produced in situ rather 
than in the more rarefied GM. This process might increase the 
species abundance, but cannot achieve the required decoupling of 
6 Li evolution from the enrichment history. 

Alternatively, shocks associated with structure forma- 
tion might provide an alternative 6 Li production channel 
llSuzuki fc Inouj |2002|) ; although potentially interesting as this 
mechanisms decouples metal enrichment (governed by SNe) and 
CR acceleration (due to structure formation shocks), the diffi- 
culties that this scenario must face are that (i) at the redshifts 
[z = 2 — 3) at which shocks are most efficient it must be still 
[Fe/H]< —3, and (ii) MPHS th at formed at ea rlier epochs should 
have vanishing 6 Li abundance jPrantzosll2006h . 

If these issues could represent insurmountable problems, 
then one has to resort to more exotic models involving either 
suitable modifications of BBN or some yet unknown production 
mechanism unrelated to cosmic SF history. 
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